N <- 100
Y1 <- rchisq(N,df=4)
hist(Y1,breaks=100,freq=F)
